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I. INTRODUCTION 



The uncertainty relations play a fundamental role not only in the foundations of quan- 
tum mechanics [1, 2] but also for the quantum description of the internal structure of d- 
dimensional physical systems [1-5] as well as for the development of quantum information 
and computation [6, 7]. The (position-momentum) uncertainty principle has attracted con- 
siderable attention since the early days of quantum mechanics [8, 9] up until now [1, 2, 10-12] 
because of its numerous scientific and technological implications. The first mathematical 
relation which expresses this principle in an exact and quantitative form is the celebrated 
Heisenberg relation [8, 9] which uses the standard deviation or its square, the variance of 
position and momentum, as measure of uncertainty; assuming (x) = (p) = for notational 
simplicity, it reads as 

{r'm > J (1) 

for d-dimensional quantum-mechanical states. 

However, this relation is not only too weak but also it is often inadequate [12-16]. In 
order to take care of these problems, various alternative formulations of the uncertainty 
principle have been proposed by use of some information-theoretic uncertainty measures 
hke the Shannon entropy [17], Renyi entropies [18-20], TsaUis entropies [21, 22], entropic 
momenta [23] and Fisher information [24-26], as recently surveyed [5, 12, 27]. 

Not so well known is the moment-based uncertainty relation developed by Angulo [28, 29] 
in 1993 which can be recast [5] under the form 

V(ae)a r3(l + gy \{he)-,V-a{l + 1)) 

valid for all (a, h) e = (0, +oo)^. This relation, which offers a more general and versatile 
formulation of the uncertainty principle (note that it reduces to the Heisenberg inequality (1) 
in the particular case a — h — 2), has not received so much attention despite the knowledge 
of the moments often completely characterize a probability density. Strictly speaking, in the 
d-dimensional case and when the characteristic function admits a Taylor expansion at any 
order, the assertion that the moments characterize a distribution is true concerning all the 

« d 

moments of the form / TT [x^^ p{-x)dxi) for all /cj e N. The assertion is no more true when 
(some of) these moments do not exist and/or dealing only with fractional moments. For 
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example, this appears for laws that are not exponentially decreasing (e.g. power law such as 
Levy noise). This is known as the Hamburger moment problem [30, chap. Ill, §8]. Finally, 
moments of various orders often describe fundamental quantities of the involved quantum 
system [5] . Other similar relationships for particular values of the parameters have also been 
published [10, 31, 32]. Note also that quantities {r"')^{p^)^ are insensitive to a stretching 
factor in the position (or equivalently in the momentum). Moreover, for specific values of 
a and/or b, the moments are linked to physical quantities (e.g. atomic Thomas- Fermi or 
Dirac exchanges [5]). Thus, it may offer a useful tool to quantify complexity for atomic or 
chemical systems that can be complementary to those proposed e.g. in [5, 33-35]. 

In this work we deal with relation (2) and improve it by use of a Rcnyi-entropy-based 
approach, in a way similar to the procedure followed by Bialynicki-Birula and Mycielski 
(BBM in short) [17] and Angulo [28, 29] used to obtain the relations (1) and (2), respectively, 
from the Shannon entropy. For this purpose we first fix notations and briefiy review the 
entropic uncertainty relations in Section II. Then, in Section III, we find a moment-based 
formulation of the uncertainty principle which extends and generalizes the relations (1) 
and (2). In Section IV we carry out a computational analysis of the new moment-based 
uncertainty relation for hydrogenic and oscillator-like systems, not only because they are 
the two main quantum prototypes in d-dimensional physics but also because their position 
and momentum moments have known analytical expressions in terms of the hyperquantum 
numbers at all orders [36] . Finally, some conclusions are given in section V. In the appendices 
we provide help to clearly discuss the proof of the moment uncertainty relation described in 
Section III. 

II. ENTROPIC UNCERTAINTY RELATIONS: A BRIEF REVIEW 

Let us denote by ^(x) and ^(p) the wavefunctions of a ci-dimensional quantum- 
mechanical system in the position and momentum spaces, respectively, so that 

$(p) = (27r)-'^/2 [ *(x) exp(-zx*p) dx 

where the units with h = I are used. The corresponding position and momentum probability 
densities will be denoted as 

p(x) = |*(x)|2 and 7(p) = |$(p)|2, 
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respectively. These two density functions are known to be completely characterized by the 
knowledge of the moments (r") and (p^) of all orders, respectively, where r = ||x|| and 
p = IIpII denote the Euclidean norms of the d-dimensional position and momentum single- 
particle operators, respectively. The position expectation value (/(r)) is defined as 

(/(r))= / /(||x||)p(x)dx, 

and similarly for the expectation value {f{p)) with the momentum density 7(p). 

For notational simplicity, we assume that x and p have zero mean, so that the variance- 
based Heisenberg uncertainty relation takes the form (1). Nowadays it is well-known that 
there exist other uncertainty relations which are much more stringent. They are based on 
information-theoretic quantities such as the Shannon and Renyi entropies and the Fisher 
information, which provide complementary measures of the position and momentum proba- 
bility spreading. Let us here recall the definition of the Renyi entropy of (real) index A > 0, 
A 7^ 1 [37, 38] 



Hxip) 



1- A 



log / (p(x))^dx 



2A 
1- A 



l0g||*||2A, 



(3) 



which represents an alternative generalized measure of uncertainty (lack of information) of 
a random variable with probability density p = Here, || • \\s denotes the norm for 

functions: ||^||p = {J^^ |^(x)|^(ix)^^*. Note that Urn Hx{p) = H{p) = - p(x) logp(x)dx 
is the Shannon entropy, that can thus be viewed as a special case of the family of Renyi 
entropies (we will write H — Hi). 

To derive an entropic formulation of the uncertainty relation, the point to start with is 
the Beckner relation that links the norm of a (wave) function ^(x) to the norm of its 
Fourier transform ^(p) where x and p are continuous in R'', d being the dimension, and s 
and q being conjugated numbers in the Holder sense: 1/s + 1/q = 1. This relation states 
that for any s G [1; 2] and q = s/{l — s): 



where 



Cs,q — 



<(<^.,.)'II*IL, 



2n 
s 



(4) 



(5) 



is the Babenko-Beckner constant. Thus, by taking the logarithm of the relation (4) with 
s — 2a and q — 2a*, one achieves the relation [18, 19] 

log(2a) log(2a* 



HM + Ha*h)>d{\og{27r) + 



2{a-l) 2{a*-l) 



(6) 



where a and a* are two real parameters related by ^ + ^ = 2, from which we define 
a* {a) — a/{2a—l). In principle a e [| ; l] but it can be seen that by symmetry (exchanging 
the roles of ^ and '^), this relation holds for any a > 1/2. When a — > 1, then a* — > 1 and 
thus the BBM relation [17] dealing with Shannon entropies is recovered 

//(p) + //(7)>d(l + log(7r)). (7) 

The entropic uncertainty relations given in (6) and (7) can be recast in the more convenient 
product form 

NMNa*il)>B{a), (8) 
with ^ ^ 

ri/ cx—l rv* a*— 1 1 

B{a) = for a ^ 1 and fi(l) = - , (9) 

using the so-called Renyi A-entropy power 

where the limiting case A ^ 1 corresponds to the Shannon entropy power N = Ni. BBM 
showed also that his primary relation (7) using Shannon entropies does imply the Heisenberg 
relation (1). To show this, it suffices to search for the maximizer of N{p) under variance 
constraint (r^) fixed, that is known to be a Gaussian of covariance matrix ^ I for which the 
entropy power is ^ [38, 39]. The same work is then done (separately) for the momentum, 
i.e. for A'"(7) subject to (p^) fixed, to finally achieve 

{r'){p')>d'N{p)N{j)>^ (11) 

and thus the Heisenberg relation. Heisenberg inequality is known to be sharp and, fortu- 
nately, nothing is lost by this way of making. Indeed, equality between the entropy and 
its maximal value is reached if and only if p is Gaussian. Furthermore, if (and only if) p 
is Gaussian, 7 is also Gaussian with the "appropriate" variance, and thus simultaneously, 
in the momentum space the maximum entropy is achieved. In other words, the sum of 
the maximum entropies corresponds to the maximum of the sum here. Simultaneously, the 
BBM inequality becomes an equality if and only if p is Gaussian, and thus the succession of 
inequalities are equalities. 

Note now that the relation (8) with Renyi entropies given above concerns only indexes a 
and a* so that 2a and 2a* are conjugated in the Holder sense: ^ + 2^ — ^- Zozor et al. 



[20] then showed that the relation (8) extends for any pair (a, /3) in such that /3 < a* (a), 
simply noting that Nx viewed as a function of A is decreasing (and after decomposing the 
allowed domain for the parameters into three regions), leading to 

NMNM>2{a,P) (12) 

where the bound is 

f l/e2 for {a,l3) e [0; 1/2]^ 

Z{a,P) = l (13) 
[ B(max(Q;, /?)) otherwise 

with B defined in Eq. (9). 

Note that on the "conjugation curve" /3 — a*{a) — a/{2a — 1), the bound is sharp and 
attained if (and only if) p is Gaussian, since it is the (only) case of equality in the Babenko- 
Beckner relation (see Lieb's paper [40]). Finally, let us also mention that Zozor et al. [20] 

showed that for /3 > a* no uncertainty principle exists, in the sense that the product of 
Renyi entropy powers is just trivially non-negative. But below the conjugation curve, it 
is not known yet neither the sharpest bound, nor the states that saturate the uncertainty 
relation. 

III. THE MOMENT-BASED UNCERTAINTY RELATION 

The uncertainty relation (2) based on the moments (r") and (p^) was obtained in [28, 
29, 41] by use of two elements: the Shannon-entropy-based BBM relation (7) and the maxi- 
mizer [41] of the Shannon entropy of the position (momentum) density subject to (s.t.) the 
constraint (r") (resp. {p^)). Let us remark that such a bound cannot be sharp. If we denote 
by ^max,a the wavcfuuctiou that gives the maximizer of N{p) s.t. (r") and by ^'max.b the 
wavefunction that maximizes A^(7) s.t. (p^), then these two functions are not linked by a 
Fourier transformation, namely ^^max.b 7^ ^max,a, except for the particular case a = b = 2. 
Or, in other words, the sum of the maximal entropies is not here the maximum of the 
sum. When deriving the Heisenberg relation from the Bialynicki-Birula relation, although 
the maximization is made separately on each Shannon entropy, it appears that the squared 
roots of the two maximizers are precisely linked by a Fourier transformation. Without en- 
tering into details here, let us consider the example of Pmax,a — argmaxi7(p) s.t. (r") that 
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is a generalized Gaussian of index a [42, 43]. Its squared root is thus a generalized Gaussian 
of index a/2 and the Fourier transform of the probability density function (pdf) is not a 
generalized Gaussian (i.e. not the maximizer of the other Shannon entropy): it is linked to 
an a-stable law of stability index a/2 [44]. 

In this Section we improve the relation (2) in a similar way but using the Renyi en- 
tropy (3), which includes the Shannon entropy as a particular case. Our procedure has the 
following steps: 

1. Start with the Renyi-entropy-based inequality (12), namely Na{p)N^{'y) > Z{a,(3), 
with the bound Z defined in Eq. (13). 

2. Search for the maximum Renyi entropy power Na{p) s.t. {r°'). This will give rise to a 
relation of the form (r")^/" > Na{p)M.{a, a), where the bound A4 has to be obtained 
in terms of a and a (see Appendices A and B). 

3. Similarly (and separately) for the momentum one will arrive at the relation {jJ^Y^'^ > 
N^{^)M{b,l3) 

4. These will lead to (r»)i(/)t > N^{p)N^{-f)M{a,a)M{b, /3) > Mia,a)M{b, /3)Z{a, /3), 
for every pair (a, 6) e M^. 

5. Finally, the best bound we can find is C(a, b) = max^^^ A^(a, a)A4{b, f3)Z{a, (3), where 
/3 < a* (a) (other restrictions come out that considerably reduce the (a,/3) domain for 
searching the maximum; see Appendix C). 

It can be shown (see Appendix C 1) that the desired maximum is on the conjugation curve 
j3 = a*{a), and then C{a, b) = max A^(a, a)M{b, a*)B{a). 

As previously mentioned, the bound must be at least the same as the case of Dehesa et al. 
[5] , since the latter corresponds to the particular situation a — P — 1 in our computations. 

The main result of the present effort is summarized here (and proved in the appendices): 

For any a > 6 > there exists an Uncertainty Principle that can be stated in the following 
way for arbitrary-order moments of the position and momentum observables in d-dimensional 
systems 

(r«)l(p'')f > C{a,b) = meixB{a)Mia,a)Mib,a*), (14) 
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where B{a) is defined in Eq. (9), a* — a/ {2a — 1), 



and the function A4 has the form 

2nel ^ 



(15) 



-rf(A-l) Y 



nB{f,i-^^-f)} \d{x-i) + ixj \d{x-i) + ix 



IX 



d(X-l) 



Mil,X) = < 



27re| 



nr{j) \ie) ' 



l-^b<^<l 



A = 1 



(16) 



27re 



d{X - 1) 



\nB{j,^)J \d{X-l) + lXj V^^(A-1) + 



IX 



IX 



d(A-l) 



, A > 1 



with fl — ^jj^ and B{x, y) the beta function. 

The case 6 > a > can be treated using the symmetry (proved in the appendix), 



Q;opt(fl)^) — argmax,B(ci;)Al(a, q;)A4(6, a*) 



satisfies 



and then 



aopt(^', a) = (Q;opt(a, 6))*, 



C{h,a) ^C{a,b). 



(17) 



(18) 



(19) 



The symmetry on a^pt allows also to conclude that Q;opt(0)Ct) — 1 and thus the optimal 
bound from our approach coincides with that of Angulo, given in (2). Unfortunately, except 
for the case a = 6, we have not been able yet to obtain an analytical expression for C{a, b). 

Figure 1 depicts the bound C{a,b) for given values of a, as a function of b, compared to 
the bound V{a,b). From the figure we see that the bound is substantially improved when 
b ^ a, especially as b departs considerably from a. 

Figure 2 depicts the optimal a — ctopt as a function of b in the same configurations as in 
Fig. 1. The curves illustrate that only for a — b, the optimal bound is obtained for ccopt = 1- 
For a ^ b, a finer study of A4 could allow to even reduce the domain D where a^pt lies. 
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5 10 5 10 5 10 5 10 5 10 

b b h b b 

FIG. 1. Bound C{a,b) (solid line) given in (14) compared to 'D{a,b) in (2) (dashed line), versus b, 
for given a = 0.1, 0.5, 1, 2 and 4 respectively, with d = 5. For each value of a, the new bound C is 
always above V; both functions coincide when b = a. 




b b b b 



FIG. 2. aopt(a, b) (solid line) given in (17), versus b, for given a = 0.1, 0.5, 1, 2 and 4 respectively, 
with d = 5. The dotted vertical line indicates b = a. Thus, left to this line, a > b > and ctopt has 
to be searched in D Eq. (15). This domain is indicated by the dashed lines. At the opposite, right 
to the vertical dotted line, b > a. Thus, symmetry Eq. (18) is used and aopt{b, a) = (aopt(a, b))* is 
seek. Since 6 > a > 0, this parameter is also in domain D Eq. (15) (where b replaces a); the dotted 
curve represents aopt{b,a) (the solid curve being aopt(a, &)) and domain D is still represented by 
the dashed lines. 



IV. APPLICATION TO CENTRAL POTENTIAL PROBLEMS 



Let us now apply and discuss the minimal uncertainty bound (14) for the two main 
prototypes of c?- dimensional physics: hydrogenic and oscillator-like systems. But before, let 
us give a brief review on eigensolutions for quantum systems in central potentials. 



A. Eigensolutions for central potentials: a brief review 



In both hydrogenic and oscillator cases, the quantum systems are described by the phys- 
ical solutions of the Schrodinger equation 

1 



-V^ + V{r) 



(20) 



where V{r) is a radial potential and where, without loss of generality, the mass is set to 
m = 1. It is well known [45] that the wavefunctions of a Hamiltonian with central potential 
can be separated out into a radial, RE,i{r), and an angular, y^ij_y{fld-i), part as 



(21) 



The position x = {xi, . . . , Xd) is given in hyperspherical coordinates as (r, 9i, 02, ... , Od-i) = 
(r, Qd-i), where naturally ||x|| — r — \J Yli=i ^1 ^ [0 j +oo) and Xi — r ^111=1 ^i^^fc) cos^j 
for 1 < i < (/ and with G [0 ; tt), i < (i — 1, 9d-i G [0 ; 27r). By convention = and the 
empty product is the unity. The angular part, common to any central potential, is given by 
the hyperspherical harmonics [45, 46] y{^i}{P'd-i)-i which are known to satisfy the eigenvalue 
equation 

ALiJ^M(^^d-i) ^l{l + d-2) J^{^}(fi<i-i), 
associated to the generalized angular momentum operator given by 



d-l 



^ (sin^i)^+i-'^ d 



1=1 



(sin 0i) 



d—i 



(n5=isin^,)' 

The angular quantum numbers {/i} = {/ii = / ,/i2 , • • • , l^d~i = ^} characterize the hy- 
perspherical harmonics, and satisfy the chain of inequalities i = /^i > /X2 > • • • > A*d-2 > 



|Afd-i| 



The radial part RE,i{r) fulfills the second-order differential equation 



1 d^ d - 1 d lil + d-2) ; 
V + Vir) 

2 dr2 2r dr 2r2 ^ ^ > 



RE,i{r)^ERE,i{r), 



which only depends on the eigenenergy E, the dimensionality d and the largest angular 
quantum number I — /ii. 

Then, the quantum-mechanical position probability density for central systems is given 

by 



(22) 
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It is worth remarking that this density function is normahzed to unity. Let us bring here 
that 

r+oo r 

/ T''~^\RE,i{r)?^r^l and / dO^.i = 1, 

Jo J[0;7r)'*-2x[0:27r) 

and that the volume element can be expressed in hyperspherical coordinates as 

dx = r^-^ dr dQd_i = r'^-^ dr ^JJ(sin O^Y'^-^ AO^j AOd-i. 
Thus, the moment (r") for the d-dimensional density Pe,{ij.}{^) has the expression 

r+oo 

^^a) = / r'^+»-i|i?^_^(r)|2dr, (23) 







which is only characterized by the position radial wavefunction Re,i{t) of the particle. 

Prom the Fourier transform of it comes out that in the momentum domain the 

wavefunction also separates under the form 

(see e.g. [36, 46, 47]) with the same hyperspherical part, and the radial part expresses from 
Re,i through the Hankel transform (e.g. [48, 49]), 

r+oo 

Me,i{p)=P^~^ r'^RE,i{r)Ji+'^-i{pr)^r (24) 

(Jjy is the Bessel function of the first kind and of order i/). Immediately, in the momentum 
space, the moment (p^) has the expression 

r+oo 

{p')= / p'+'-'\Me,i{p)\' dp, (25) 
Jo 

which is only characterized by the momentum radial wavefunction Me,i{p) of the particle. 

These expressions have allowed to find numerous information-theoretic properties [24, 
25, 36, 45, 50, 51] of general central potentials, particularly the Heisenberg [25] and Fisher- 
information [24, 25] uncertainty relations, as recently reviewed [5]. 



B. Application to d-dimensional hydrogenic systems 



Let us now examine the accuracy of the moments-based uncertainty relations (14) for 
the main prototype of d-dimensional systems, namely the hydrogenic atom. This system 
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has been recently investigated in Ref. [36] in full detail from the information theory point 
of view. In this case, the potential has the form V{r) — (without loss of generahty, the 
atomic number Z is taken to be 1) and the energies are 

where rj denotes the grand principal quantum number. The radial part of the eigenfunctions 
are thus completely calculable [36, 50, 51]. The radial wavefunction in position domain 
expresses as 



where L = I + / = 0, ... ,n — 1 is the grand orbital quantum number, f = ^ is a 
reduced (dimensionless) position and £^ are the Laguerre polynomials. As it is shown in 
Refs. [36, 50, 51], after the Hankel transform (24), the radial wavefunction in momentum 
domain expresses as 



^^'^^^^ ^ ' V 2vrr(r, + L + l) + ^ ^ (iT^ ^-^"^ [iTW ) 

where p = Typ is a reduced (dimensionless) momentum and are the Gegenbauer polyno- 
mials. Prom these expressions together with (23) and (25), it is shown [36] that the position 
and momentum moments of arbitrary orders, corresponding to a given eigenstate character- 
ized by an energy E and an angular quantum number / (or equivalently by r] and L) have 
the expressions: 

^""^ = '^°2"+irf2L 3i^2(-ry + L + 1, -a - 1, a + 2 ; 2L + 2, 1 ; 1) 



and 



(27) 



^ 4r(ry + L + l)r(L+^) r(L + ^) 
^ ^ v'-' r{rj - L) r2 (L + 1) r(2L + 4) 

/ 6+3 5-6 3 5 

5F4fL-?7+l, L+77+1, L+1, L+ — , L+— ; 2L+2, L+-, L+2, L+-; 1 

for 6 < 2L-I-5, where pFq are the generahzed hypergeometric functions (see e.g. [52, 2.19.14- 
15] and reflective properties of hypergeometric functions). Note that the momentum wave 
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function is not exponentially decreasing. The direct consequence is that not all moments 
exist in the momentum domain, what is reflected in the restriction for the values of b. 

Thus, the uncertainty product (r")a(p*)§ can be computed and therefore studied analyt- 
ically for hydrogenic systems in c?- dimensions. As an illustration. Fig. 3 depicts the product 
(r")^(p^)f computed from (27), for (a, 6) = (1,2) and Fig. 4 plots the case (a, 6) = (1,4) 
(both for d = 3) when the system is in the state {E,l), together with the corresponding 
bound C{a,b) given in (14)-(16). 

j ° ° o 

o o ° 

° o ° 

_ ° o ° 

'^^ 1 o 

= ° 

o 

□ □□□□□□□□□□□□□□ 



n=l n=2 n=2 n=3 n=3 n=3 n=4 n=4 n=4 n=4 n=5 n=5 n=5 n=5 n=5 
1=0 1=0 1=1 1=0 1=1 1=2 1=0 1=1 1=2 1=3 1=0 1=1 1=2 1=3 1=4 

FIG. 3. Product (r)^(p^), i.e. (a, 6) = (1,2) in eqs. (27) (circles) for the lowest energy states and 
lower bound C(l, 2), eqs. (14)-(16), of this product (squares), for 3- dimensional hydrogenic systems 
(d = 3). 



\ \ I \ \ \ \ I \ \ \ \ I \ \ 
n=l n=2 n=2 n=3 n=3 n=3 n=4 n=4 n=4 n=4 n=5 n=5 n=5 n=5 n=5 
1=0 1=0 1=1 1=0 1=1 1=2 1=0 1=1 1=2 1=3 1=0 1=1 1=2 1=3 1=4 

FIG. 4. Product {r)^{p^)2, i.e. (a, 6) = (1,4) in eqs. (27) (circles) for the lowest energy states and 
lower bound C(l,4), eqs. (14)-(16), of this product (squares), for 3-dimensional hydrogenic systems 
(d = 3). 

We can see from both flgures that, although not sharp, the bound C{a, b) is close to the 
product (r°)o(p^)i for the ground state (n = 1 and / = 0). However, when n increases, the 
discrepancy from the bound increases (and decreases with / for flxed n). The same behavior 
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occurs for other pairs (a, b) whatever the dimensionahty d. Since hydrogenic systems belong 
to the family of radial potential systems, this suggests that a refinement can be found 
in the context of radial systems as already done for the usual variance-based Heisenberg 
inequality, and for Fisher information-based versions [24, 25]. To give a further illustration. 
Fig. 5 depicts {r"')'^{p^)i as a function of 6, for fixed a, and for the ground state {n = 1, 
/ = 0) in 3 dimensions. In all the cases shown, we observe the existence of a value of h that 
is "optimum" in the sense that the uncertainty product is close to the bound proposed here, 
corresponding to a situation of low generalized uncertainty. As h increases (up to 2L -|- 5 = 5 
for the ground state in 3 dimensions), the uncertainty departs from our bound. Finally, one 
observes for the tested values of a that the lower bound has a concave behavior versus 6, 
while the product (r'^)a(p*)^ exhibits a convex behavior. This suggests the existence of an 
optimal value of h (function of a) in terms of low discrepancy from the bound. 
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FIG. 5. Product (r")a(p^)b (solid lines) in the ground state (n = 1, / = 0), for fixed a = 0.1, 
a = 0.5, a = 1, a = 2 and a = 4 respectively, and lower bound (dashed lines) for 3-dimensional 
hydrogenic systems (d = 3). 



C. Application to d-dimensional oscillator-like systems 

Let us consider now a potential of the form V{r) = (without loss of generality, the 
product mass squared pulsation is taken as mu'^ = !)• In this case, the energies are 

E = 2n + l + ^, 72 = 0,1,... and / = 0, 1,... 

and the radial parts of the wavefunctions are again known [53]. They express as 
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and Mej{p) = Re,i{p)- Comparing (28) with (26), after a change of variables r = r^, one 
can easily show from (27) that the statistical moments write down as 

(r^) = l-A 3F2 -n, , - + 1 : L + -, 1 : 1 

\ / Y{l + I) \ ' 2'2 ' 2' ' 

and (29) 

r(/ + ^) / b b d 

(/)-^^(^3F.(-n,--,- + l;L+-,l;l 

(see also Ref. [54] for special cases). 

Figure 6 describes the moments product {r"-)^{p^)T> using (29), for (a, 6) = (1,2) and 
Fig. 7 exhibits the case (a, 6) = (1,4) (both for d = 3), together with the corresponding 
bound C{a,b) given in (14)-(16). 

= 00° 

o ° 



^10 



o 

o 



o o 

o 



\ I \ \ \ \ \ \ \ I \ \ \ \ \ I 
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FIG. 6. Product (r)^(p^), i.e. (a, 6) = (1,2) in eqs. (27) (circles) for the lowest energy states and 
lower bound C(l, 2), eqs. (14)-(16), of this product (squares) for 3-dimensional harmonic oscillators 
(d = 3). 

We can see from these figures also that even if not sharp, the bound C{a, b) is very close 
to the product (r")a(p^)5, for the ground state (n = and / = 0). The global behavior is 
similar to what happens for the hydrogenic systems: there is a discrepancy from the bound 
as n increases. But here, the discrepancy increases also with / when n is fixed. The same 
behavior occurs for other pairs (a, b) and whatever the dimension d. In fact, when observing 
more finely {r)'^{p'^) and {r)'^{p^)^, it appears that these products essentially depend on the 
energy level, i.e the values of these products for a fixed value of 2n + Z are very close (see 
e.g. n = 0, / = 2 or n = 1, / = 0). This was true also for the hydrogenic systems, but it is 
strongly more pronounced for the harmonic oscillator. All these observations reinforce our 
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FIG. 7. Product {r)^{p^)^, i.e. (a, 6) = (1,4) in eqs. (27) (circles) for the lowest energy states and 
lower bound C(l, 4), eqs. (14)-(16), of this product (squares) for 3-dimensional harmonic oscillators 
(d = 3). 

"conjecture" that refinement can be found in the context of radial systems, for moments' 
orders other than a = 6 = 2, at least in terms of energy levels. 

A further illustration is given by Fig. 8 where (r")a versus h is depicted, for different 
fixed values of a, in the case of the ground state (n = 0, / = 0). 
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FIG. 8. Product (r")o(p'')fc (solid lines) in the ground state (n = 0, / = 0), for fixed a = 0.1, 
a = 0.5,a = l,a = 2 and a = 4 respectively, and lower bound (dashed lines) for 3-dimensional 
oscillators (d = 3). 



Globally the behavior of the moments' product compared to the bound observed here is 
similar to that of the hydrogenic systems. However, the discrepancy from the bound is less 
pronounced for the harmonic oscillator (in the ground state) than for the hydrogen systems. 
Note that the bound is achieved in the case when a = h = 2. This case corresponds 
to the classical variance-based Heisenberg inequality. Moreover, the ground state of the 
oscillator leads to the Gaussian pdf p (and 7): in this case the variance-based Heisenberg 
inequality is saturated. One can again observe the convexity of the product (r'^)a(p*)f (in 
fact almost linear): together with the observed concavity of the lower bound, this reinforce 
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our conjecture on the existence of an optimal value of 6(a) in terms of low discrepancy from 
the bound. This remains to be studied more systematically and more deeply. 



V. DISCUSSION AND CONCLUSIONS 



In this paper we have proposed an improved version of the moment-based mathematical 
formulation of the position-momentum uncertainty principle for quantum systems that gen- 
eralizes the seminal variance-based formulation of Heisenberg. The main result of this contri- 
bution is formalized in Eq. (14) together with eqs. (15), (16) and (9): (r")a(p'')l > C{a,b). 
In contrast to the entropic uncertainty relations (like Eq. (12)), the new formulation is 
based on spreading measures which describe physical observables. Our present approach 
suffers, however, from the fact that the lower bound C{a,b) found here for the product of 
the position and momentum moments for arbitrary a and h is not sharp. To tackle this issue 
a variational approach may be envisaged, although it is a difficult task. Another alternative 
might be to employ appropriate Sobolev-like inequalities, as done for entropic formulations 
(see e.g. Ref. [17-19, 21]). 

The new moment-based uncertainty relation is physico-computationally analyzed in some 
d-dimensional quantum systems. Precisely, the bound of the moment-based uncertainty re- 
lation is compared to the product of the moments for hydrogenic and oscillator-like systems. 
In both cases analytic expressions of the moments exist in terms of hypergeometric functions 
(eqs. (27) and (29) respectively). Our results suggest that the improvement of this relation 
for general central potentials seems possible whatever the orders a and h of the moments, at 
least in terms of energy levels. Such an improvement exists in the variance-based context 
a = b = 2 [24, 25], but for moments of arbitrary order this issue is a fully open problem 
which deserves to be variationally solved for both fundamental and applied reasons. This 
suggests also that the product (r") « (j)^) 5 can be envisaged as a useful tool to quantify the 
complexity and organization of various physical systems. However, the properties of such a 
complexity measure should be analyzed in more detail. 
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Appendix A: Evaluation of the maximizers of the Renyi entropy power under mo- 
ment constraint 

In section III (steps 2 and 3) we established the necessity of searching for the maximum 
of the Renyi entropy power Na{p) subjected to given moment (r") and Nj^l'-f) s.t. (p^). This 
variational problem has been tackled and partially solved by Dehesa et al. [55]. Similarly 
to what is done for variance constraint, the problem is to maximize the frequency entropic 
moment of order A > 0, an increasing function of the entropy power, / /(x)'^dx s.t. 

/ /Wd. . : and / - (.'), w:. , > a., w.e. - r lie E.Mean 

norm of x. Note that we work here with the variables x and r, but the results obtained 
will be valid in the momentum domain, changing x to p and r to p. Then, we have to 
maximize / (/(x)^ — ///(x) — i/||x||'/(x)) dx, where // and v are the Lagrange factors. 

JRd 

Prom the corresponding Euler-Lagrange equation, one obtains that / must be of the form 

/(x) = ^^^±^dWLj ^ ^ ^ where (t/)+ = max(|/,0). With integrability arguments (/ must be a 
pdf, and thus positive and integrable), /i > and u must have the sign of 1 — A, and thus 
the pdf that maximizes the entropy power Nx s.t. (r') can be recast under the form 

/,,(x) = C(l-(A-l)(||x||/5)')^. (Al) 

This pdf is sometimes called generalized Gaussian [56, 57], but this terminology is not 
adequate. Indeed, when A — )> 1, this pdf tends to /ij(x) = Cexp(— ||x/5||^) that is also 
sometimes named generalized Gaussian (or also Kotz-type) [43, 58, 59] (and also sometimes 
as stretched exponential or power exponential [59, 60]). Purthermore, when I — 2 one 
can recognize in (Al) the well known g-Gaussian (also known as Student-t or Student-r 
depending on the sign of 1 — A), where g = 2 — A and thus the generalization (Al) is known 
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under the terminology of stretched q-exponential [55, 61] or even generalized q-Gaussian of 
parameter q — 2 — X and (stretching) parameter I. 

Constants C and 5 are to be determined so that the constraints are satisfied. The 
normahzation constraint reads 



C I (l-(A-l)(|M|/5)')r dx 

Cn r"^-^ {1- {X-l){r/Sy)^-' dr 
Jo 

^ I i'^/'-^ll-signlA-ljt)!- dt 



where the second hne comes from [62, Eq. 4.642], with Vt = ^(1/2) ^^^^ ^hc surface of 
the (i-dimensional unit sphere and where the third hne comes from the change of variable 
r = ^{t/\i — The integral term, that we will denote Bi{l,X), expresses via the beta 

function B{x,y) = r{x)r{y)/r{x + y), from [62, Eq. 8.380-1 & 8.380-3], one finally obtains 

^~ l\X-l\'^/^ 

where 



Bi{l,X) (A2) 



Biil,X) = { (A3) 

Indeed, the integral converges provided that A > 1 — l/d. 
In the same vein, the power moment constraint writes 

{r^)^C I ||x||'(l-(A-l)(||x||/5)0f^dx 

= Cn r^+'^-^ (1 - (A - l)(r/5)')+"' dr 
Jo 

- ^|A-i|./m I (1 - «ig-(^ - 1) ^y^-' d^ 
where the integral term, denoted here 3^(1, X), expresses via the beta function from [62, 
Eq. 8.380-1 & 8.380-3], leading to 

fi nxi+d 

(^') - JlX^TW^l^rnilA) (A4) 

where 

, , + x^) if A>1 

B^{l,X)={ (A5) 
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Note that the existence of the latter integral implies a stronger restriction to A than the one 
coming from the normalization, that is we require now: 



In both constraints, the case A = 1 can be recovered by letting A — > 1+ or A — 1~ (from 
[63, 6.1.47] or [62, 8.328-1], limy^^ B{x,y)x^ = r(x)): it is not needed to treat this case 
separately. 



Appendix B: Maximal entropy power Nx and bound for the moment (r ) 

Following the procedure proposed in section III, we discuss here the bounds for the 
moments (r") and (p''). Prom (Al), the maximal A-norm of fx,i{x to the power A takes the 
form 



ll/A,dlA = C'^ / (l-(A-l)(||x||/5y)^-dx 

P + OO 

^C^Q r'^-^l-{X-l){r/Sy)'-' dr 
Jo 

/ f^/^-' (1 - sign(A - 1) t)^-^ dt. 
Jo 



/lA- IK/' 



Then, from [62, Eq. 8.380-1 & 8.380-3] we obtain 



fix QXd 

ll/A,dl^- .^_ '/^ ^/.(U), (Bl) 



where we have defined 



, , , + if A>1 

Bh{l,X)-{ (B2) 
i^-f) if 1-^<A<1 



20 



that adds no new restriction on A. Thus, the maximal value of the Renyi entropy power is 
Nx{h,i) 



2 

1 / — \ 



2ne 



A 



27re I yilX-ll"^/^ 
27re I V ^1 



from (A2) and where the arguments of Bi and are omitted for simplicity. Taking the 
ratio ''^a/i+i from (A2) and (A4), one obtains 

d/l 



2 



that gives 



One can simplify a little bit this expression by considering the parameter 

= /.(A) = ^ (B3) 

that governs the maximal entropy power, with > 1 or < —d/l. Noting that 

Bi d + li^ Bh III , . 

— = sign(;.)— and ^ = ^ (B4) 



2 



SO that 



27re I I \siga(ijL) d J \ l/j, 

We finally obtain that the Renyi entropy power of any pdf p, Nx{p) — ^ ^IIpIIa"^ J > s.t. 
fixed (r'), is bounded from above by the maximum value Nx{f\,i). Therefore we can write 



{r'y/'>N,ip)M{l,X) (B6) 
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where function Ai expresses as 
27re ' ^ 



nB{f,ij)J \d + liij \d + liJ, 



2(m-1) 



M{l,X) = < 



27re 



27re 



d\^ 



2 ^ ^ ''(M-l) 

d \' f //i ^ 



if A > 1 



if A = l 



with 



A 



A- 1 



(B7) 



(B8) 



and where J\4{1, 1) = hm^^i A4{1, A) from the first and/or second expression of A4 and [63, 
6.1.41]. 



Appendix C: Generalized Heisenberg-like uncertainty relation 

Using (B6) apphed to r with I — a and X — a and apphed to p with I — b and \ — j3 
respectively, and using (12), we achieve the relation estabhshed in point 4 of section III 

(Ot (jl')l > Z{a,P)M(a,a)M(b,P) (CI) 

for all a, 6 > 0, q; > /3 > /3 < and with the bounds Z and B given in eqs. (13) 
and (B7). 



1. The maximal bound is on the conjugation curve f3 = a* 

We will now show that the pair {a,f3) that maximizes Z{a, l3)Ai{a,a)Ai{b, (3) is on the 
conjugation curve, namely for 13 = a* = a /{2a — 1), for any values of a and b (under the 
existence condition for A4). 



a. Function M(l, A) is increasing with A 

Let us first consider the derivative of A4{1, A) versus A. 
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For A > 1, i.e. = ^ > 1, 



logA^ = — --logr(/x) + -logT U+ 7 - 7log(rf + /;u) + — ^log ' ^ 



dfi \ d d \ I J I d yd + l/j, 

^ '-V'(At)+V'(/x+y ) +-^-^-- + log' 



d\ \ IJ d + l/i II \d + lii 

where ip{x) = ^logr(a;) is the digamma function. 
Similarly, for A e (l — l), i.e. < —d/l, 



d 5/2 / d\ 2 

logM = ^ -^logT 1 -/i- - + -logr(l ^/i) 



dfjL dfjL \ d \ I J d 

- - log(-o? - III) + — ^ log ' 



I d \d + lfi 

^(-^(l-^)+^(l-^-^)+^-i + log(^) 



the last simplification coming from [62, Eq. 8.365-1]. 

To summarize, noting that d/i/dX = — 1/(A — 1)^, 

£logM{l,\) = ]_ ) + ^ - log/. 
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l^logM^A) = ^(A - 1)2 (^(-^) - ^"g(-^) 



( - 7 ) + log ( -/^ - 7 ) ) if 1 - < A < 1 



(C2) 

Taking the limit A — > 1"*" from the first line, or A — > 1~ from the second line, and M. being 
continuous in A = 1, we achieve ^ log A^(Z, A) |^^^ = j. 

Let us consider now the terms in the parentheses in the right-hand side of the first line 
in Eq. (C2). They can be written as g{^) — + d/l) with, 

^(/i) = V(/^) + --log//. (C3) 
A* 
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Then, from [63, 6.4.1] 



r+oo J. r+oo r+oo 

= / -e-^^dt- / e-^^dt- / te-i'^dt 

Jo 1 - e-* Jo Jo 

Jo 



-+00 _i I I ^ -t 



1 - e-* 

Now, it is easy to show that — 1 + e~* + ie~* < for i > that permits to conclude that 
g' < and thus that g is decreasing. As a conclusion, g^j,) — g^ji + d/l) > and thus 
^ log > 0: is increasing in (1 ; +00). 

Similarly, the terms in parentheses in the right-hand side of second line {/i < —d/l < 

here), writes /i(/i) — h{ix + d/l) with, 

h{^J) = '0(-/^) - iog(-//) 

and give from [63, 6.4.1] 

A* 

"+00 ^ /'+00 

ei'^dt 



r+oo J. r- 



/ — —e^^dt. 

Jo 



1 - e-* 



Then, it is easy to show that 1 — i — e"* < for i > that permits to conclude that h' < 
and thus that h is decreasing. As a conclusion, /i(/x) — h{n + d/l) > and thus also for 
^ (1 ~ /ipd j 1) 'w^^ have log > 0: is increasing. 



b. B{X) increases with AG [1/2 ; 1] and decreases with A > 1 



Prom (9) and A* = A/ (2 A — 1), the derivative of B{X) writes 



d , d /logA logA* 
logS(A) = — - + 



dX ° ' ' dX \X-l X* -1^ 

9 / log A \ d / log A* A dX* 



dX\X-lJ dX* \X*-lJ dX 

1 logA \ / 1 log A* 



^2 



A(A-l) (A-l)V \X*{X*-1) (A* - 1) V (2A - 1) 

that is 

^ log B{X) = — ^ - I - log(2A - 1) ) . (C4) 



dX ° ' ' (A - 1)H A 
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A short study of the right hand side, show that this quantity is positive if A e [1/2 ; 1] and 
negative if A > 1: B increases with A in [1/2 ; 1] and then decreases for larger values of A. 



c. Domain where the maximal bound has to be searched 



Remind that starting from (CI), namely (r")a {p^)b > Z(a, (3)Ai{a, a)Ai{b, (3), the best 
bound is then so that Z{a, l3)A4{a,a)Ai{b, f3) is maximized as a function of a and /3. Let 
us now consider the following sets in the {a, j3) plane: 

' = {{a,/3) e M^l a > < a*} 

Df, = {{a,l3) eRl\P>l,a<P*} 

Sa - {(«,/?)€ [0; l]2|/3<a} 

= {K/3)e[0; l]></3} 

where a* — and /3* = These sets are represented in Fig. 9. 

/3 1 




FIG. 9. Sets Da, Dp, Sa and Sf^ in the plane (a,/?). The solid curve represents the pairs of 
conjugated parameters, i.e. /3 = a*. The dotted arrows indicate that Z{a, P)M.{a,a)M.{b, 13) 
increases when the pair (a, (3) moves along their directions, in the sets where they are plotted. 

In order to study the best bound, we consider each subset: 

• We first consider domain and fix a. From Eq. (13), the bound is then 

Z{a, p)M{a, a)M{b, — B{a)M{a, a)M{b, P) and from the previous study of M 
we can know that it increases with f3. Thus, the bound is maximum precisely on the 
conjugation curve (3 — a*. 



25 



• By symmetry, in domain Djs and fixing one shows again that the bound is maximal 
on the conjugation curve a — (3*. 

• In the domain Si we discuss the following cases: 

— The maximum bound must be achieved on the line segment a = ^. Indeed, in Sa, 

the bound is given by B{a)M{a,a)M{b, f3) if a > 1/2 and M{a, a)M{b, f3) /e^ 
otherwise. Again, fixing a, the bound is increasing with /3 and thus is maximum 
foT (3 — a. This remains valid, by symmetry, in Sjs, and thus in all Si. 

— For a < 1/2, on the line segment a — /3 the bound is M.{a,a)M.{b,a)/e'^ and 
thus increases with a: it is maximum for a = 1/2. 

— For a G (1/2 ; 1], the bound expresses B{a)Ai{a,a)Jli{b,a) and tends to 
M.{a,a)M.{b,a)/e'^ for a 1/2. B being an increasing function in [1/2 ; 1] 
and since A4 is increasing, the bound is then maximum for a = 1. 

As a conclusion, on 5*1 the maximum bound is achieved when a = P = 1 that is again 
on the conjugation curve. 



2. Maximal bound and properties 

The best bound of the generalized Heisenberg relation one can achieve by our approach 
is then 

C{a,b)= max B(a)M{a.o)M(b,a*) (C5) 

aeD{a,b) 

where the domain of search D is ruled by the restriction on the domain of existence of Ai. 
We will come back later on this domain. 



a. Symmetries 

Let us denote by aopt{o,, b) the index that leads to C(a, 6), i.e. 

otopt{0',b) — argmaxB(Q;)A4(a, q;)A4(6, a*) (C6) 

a 

Noticing that B{a) = B{a*), one immediately observes from (C6) that 

C{b,a) = C{a,b) 
aopt{b,a) = {aopt{a,b)y. 
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(C7) 



Thus, without loss of generahty, one can restrict the study to the case with a > b. 



b. Reduced domain of search 
Consider the situation where a> b. 

If « > 1, then a* < 1. We will show that the bound B{a)A4{a,a)A4{b,a*) decreases 
with a; thus the maximum must satisfy a < 1. 

We have already seen that B{a) decreases when a > 1. Consider then the part 
M.{a,a)M.{b,a*). Remembering that a* — 2^^, one has ^ = ~J2^^^- Moreover, 
one has (-r^ = j^Ejf and ii{a*) = = --^ = -A*(q;) from Eq. (B3). Then, 

from (C2), 

d d da* d 

— log {M{a,a)M{b,a*)) = — log M{a, a) + ——log M{b, a*) 

2 /I 1 

= M TTT ~ '^il^ + ^«) + ^(/^ ~ ^b) + log(/x + da) - log(/x - db) 

where /i stands for ii{a), da — d/a and d^ — d/b. The goal is then to show the negativity of 

4, 4) = ^—r - i^il^ + da) + log()U + da) + ipiiJ, - db) - log(^ - db) (C8) 

II 11 + da 

keeping in mind that da < db- To this end, we can view this function in terms of da for 

instance, and thus the sense of variation of k{fi, da, db) = —g{fi + da)+g{fJ^ — db) + l/li — l/{n — 

db) is the same than the sense of variation of —g{fj. + da) = —'ipifJ' + da) — + log(/i + da) 

introduced Eq. (C3), versus da- We have shown that function g is decreasing and thus k is 

increasing with da- Since, da < db, to show that k is negative, it is then sufficient to show 

that k{ii,db,db) is negative. Prom [63, 6.4.1] 

dk{n,db,db) 1 1 1 

adb fJ' + db [11 + db)^ fJ'- db 



hoc 



1 + t ^—] e-""' + ( 1 ^— ] 6+""' 

1 — e~* / V 1 — e" 



/ [(1 - e-* - te-') e-""' + {l - t - e"*) e+'^'*] 
Jo 1 ~ 



dt- 



Now, it is quite easy to show that the term in square brackets is decreasing with db since 
the derivative in db is negative (the factors of e^'^^^ are negative). For = 0, it is no more 
hard to show that the square bracket is negative, that permits to conclude that for any db 
the square bracket term is negative: k{ii, db, db) decreases with db- 
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At last k{fj,, 0, 0) = and thus k{ii, db, db) < 0, implying that for any da < di, one has 
da, db) < 0. 

As claimed, -^\og {M {a, a) M{b, a*)) < for a > 1. Together with the decrease of B 
when a > 1 wc conclude that the maximum of B{a)Ai{a,a)Ai(b,a*) is attained for a < 1. 

To finish to determine the domain of search for the maximal bound, when a > b, implies 
that a < 1 as just found, a > 1/2 since it must by on the conjugation curve, and from (C2), 



In summary, for a > b, 



C(a,6)= max B{a)M{a,a)M{b,a*) (C9) 



where B and M. are respectively given by (9) and (B7). 
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